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APPLICATION OF PHASE -CHANGE TECHNIQUE TO THIN SECTIONS 
. WITH HEATING ON BOTH SURFACES 

By James L. Hunt, Joan I. Pitts, and Christine B. Richie 
Langley Research Center 

SUMMARY 


A numerical digital computer, program has been developed to calculate the heat- 
transfer coefficients for both sides of a one -dimensional finite slab subject to the bound- 
ary conditions ascribed to the phase-change coating heat -transfer technique. In a typical 
tunnel test situation where a thin wing was exposed to heating on both sides, die data 
reduction procedures for a semi-infinite slab gave heat-transfer coefficients as much as 
375 percent too high on the side with the lowest heating. The results from the one- 
dimensional finite-slab procedure are presented in the form of correction factors to the 
solution for a semi-infinite slab in terms of parameters normally used with the phase- 
change heat -transfer technique. These correlations are not restricted to slab thickness 
or thermophysical properties and are easily used to obtain accurate data on thin model 
sections. 


INTRODUCTION 


Existing solutions (refr 1) for obtaining quantitative heat-transfer data from the- 


phase -change coating technique are restricted by the assumption of a one -dimensional 
semi-infinite slab. In practice this assumption requires that the depth of heat penetra- 
tion into the model be small compared with the wall thickness. The depth of heat pene- 
tration has been shown to depend primarily upon the thermal properties of the model • 
material and the test time. A simple empirical equation given in reference 1 as 


c 0.21 2 

d ,1 a 


(1) 

/ . 


where l is the thickness of the slab and a the thermal diffusivity, determines the 
maximum test time for which the assumption of a semi -infinite slab can be applied with 
negligible error. However, models of many configurations tested in hypersonic wind 
tunnels have very thin wings and tails so that practical test times exceed the allowable 
time given by equation (1). In addition, for many tests, heat enters at both surfaces of 
such thin wings. At the present time no solutions are available for determining the 


magnitude of the errors incurred when the semi-infinite-slab assumption is violated. 
Therefore, a computer program was developed to calculate the heat -transfer coefficients 
for both sides of a one -dimensional finite slab subject to the boundary conditions ascribed 
to the phase-change coating technique (ref. 1). The results are presented in the form of 
correction factors to solutions for a semi-infinite slab in terms of parameters normally 
used with the technique. These correction factors are not restricted to slab thickness or 
thermophysical properties and can be easily used to obtain accurate data on thin model 
sections. 
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SYMBOLS 

area 

specific heat 

aerodynamic heat -transfer coefficient 

heat-transfer coefficient to stagnation point of reference sphere 
thermal conductivity 


N 

thickness of slab, V w 


l 

j=l 


'i 


free-stream Mach number 

number of blocks upon which heat balance is performed 
finite-slab heat-transfer parameter (Nusselt number), ^ 
heating rate 


q 


heating rate per unit area 


R free-stream Reynolds number 

T temperature 


T 


temperature parameter, 
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time 

thermal -interference diffusion time for thickness of slab 

thermal -interference diffusion time for one-half thickness of slab 

delta time or computing interval 

volume of block 

thickness of block 

coordinates 

chord length 

thermal diffusivity, — 

heat-transfer parameter for semi-infinite slab, 

emissivity 

angle of attack 

density " ■ - - - ; . - - - - 

Stefan-Boltzmann constant 

/yf 

time parameter for finite slab, 

7^ 


Subscripts: 

aw adiabatic or recovery conditions 

BG background 

c convection surface 


i 


initial conditions 



j variable integer (either 1 or 2) which refers to surface of slab at which heat 

is entering 

m variable integer which refers to block upon which heat balance is being 

performed 

n variable integer which refers to any other block that affects heat balance 

for block m 

pc phase change 

sis semi-infinite slab 

w wall 

1 front side of finite slab 

2 back side of finite slab 

A prime indicates the value in the previous time step. 

ANALYTICAL METHOD FOR FINITE SLAB 

A computer code (appendix A) which calculates the heat -transfer coefficients for 
both sides of a one-dimensional finite slab subject to the boundary condition ascribed to 
the phase-change coating was developed in the following manner. The finite slab of 
thickness l (fig. 1) with heat exchange at both surfaces is divided into N blocks. The 
block thicknesses are graduated inversely with the distribution of the temperature gradi- 
ent expected in the slab. These thicknesses are given in appendix A as a percentage of 
the total thickness of the slab l for two values of N (20 and 30). These two thickness 
distributions are built into the computer code and selected according to the input value 
of N. A heat balance is performed on each block. In the heat-stored term, each block 
is assumed to be at a uniform temperature throughout its volume. In the heat-conduction 
term, a linear temperature variation between centers of blocks is assumed, with heat 
passing through an area equal to the touching surfaces. 

Convection blocks (surfaces 1 and 2): 

^convection = ^j A c,j = h j A c,j( T aw,j " T c,j) 
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^radiation CT£ j A c,j[( T c,j) “ ( T Bg)J 

a v ~ 

^stored ” ^m c m v m At 

q = 2k A Tm ~ Tn 

^conducted m,n m,n w m + w n 

^convection ~ ^radiation ” ^stored ^conducted “ ® 

The surface temperatures T • are assumed to be equal to the temperature of their 

J 

respective surface blocks T £ m , which are kept extremely thin. 

Conduction blocks (internal): 

T - T 1 

n v m m 

^stored _ ^m c m v m At 

T - T 

6 = 2k A m n 

^conducted m,n m,n w + w 

’ ’ m T n 

^conducted-in ” ^stored ” ^conducted-out " ® 

“ The resultingsysterrrof linear equations is solved simultaneously at each succes- — 
sive time step with a subroutine (ref. 2) which solves the matrix equation AZ = B 
(Z denotes the unknown variables taken as temperatures of blocks here, A denotes a 
~ square coefficient matrix, and B denotes a vector of constants). The values of specific ~ 
heat and thermal conductivity are updated at the new temperatures obtained for each 
block after each time step. This heat-balance procedure is implicit and is therefore not 
bound by a stability criterion. The resulting computer program was checked against 
exact solutions of finite-slab problems (ref. 3) to determine the effect of the number of 
blocks, the distribution of the thickness of the blocks, and the computational time 
interval (At) on the accuracy of the solution for various heating rates, slab thicknesses, 
and thermophysical properties. (See appendix A.) 

As described thus far, this program calculates the temperature -time profile in each 
block for a prescribed heat flux q^ = hj(T aw ^ - T c j) and = h2( T aw 2 ~ ^c 2 )’ 

This procedure is now incorporated into an iterative cycle which enables the determina- 
tion of the heat-transfer coefficients hj and h 2 - This iterative cycle includes (1) a 
method to approximate the heat-transfer coefficient at the two surfaces to be used as 
starting points in the cycle and (2) a method of searching for the actual heat-transfer 
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coefficients - defined here as that combination which produces surface -temperature - 
time profiles that include the corresponding input points (Tp C j, tp C> j) and (Tp C 2, ^-pc, 2 ) 

(T pc i not necessarily equal to Tp C 2 and tp C 1 not necessarily equal to tp C> 2)- 

The solution to the one -dimensional heat-conduction equation for a semi-infinite 
slab with a step input in h is used to obtain the heat -transfer coefficient at both sur- 
faces for the starting point in the iterative cycle. The solution is given in appendix B. 

The approximation for a semi-infinite slab will always give a higher value for the 
heat-transfer coefficient at either surface of the finite slab than that which actually 
occurs, provided that heat is entering the slab at both surfaces (no negative heating 
rates). This is because the larger effective mass of the semi-infinite slab requires 
more heat input to reach a given surface temperature than does the finite slab. When 
the heat -transfer coefficient for the semi-infinite slab is imposed on the finite slab, the 
resulting surface temperature is higher than that which actually occurs on the surface of 
the finite slab (see sketch). As shown in the sketch, the heat-transfer coefficient is 



Surface I 


Surface 2 


Example temperature -time profile at both 
surfaces of finite slab 


higher for surface 1 than for surface 2 . The task is to find the combination of hj 
and h2 which corresponds to the respective set of surface temperature -time profiles 
(labeled "actual" in the sketch) that include the respective input points (Tp C> i, tp C ^ 
and (Tp C> 2> tpc, 2) • 
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A Newton-Raphson procedure (ref. 4) is used to determine the respective amounts 
by which both surface heat-transfer coefficients are changed at the beginning of each suc- 
cessive iterative cycle in order to reach the actual set of heat-transfer coefficients for 
the finite slab in the fewest steps. This procedure is employed at the constant input time 
points t = tp C i and t = tp C 2 and proceeds until a combination of heat-transfer coeffi- 
cients is found (hi and I 12 ) that produces corresponding surface temperatures within a 
given AT of the respective input temperatures Tp C 1 and Tp C 2 - 

RESULTS AND DISCUSSION 


Experimental Test Situation 


As an example of the magnitude of errors which may occur by violating the semi- 
infinite-slab assumption in practical test situations, phase-change data (ref. 5) obtained 
in the Langley Mach 8 variable -density hypersonic tunnel (on a 0.005 -scale Sty cast model 
of a NASA Lyndon B. Johnson Space Center space shuttle launch configuration at zero 
degrees angle of attack) were examined. The wings of the orbiter of this configuration 
were very thin (midsection thickness on the order of 0.2 cm). A chordwise profile of the 
wing on the orbiter at 69.3 percent semispan is shown in figure 2. The time at which the 

phase change (339 K) occurred along the top and bottom surfaces at this span position and 

\ 2 \ 


the thermal -interference diffusion time 
-of-the wing profile are included in “figure 27 


_ 0.2 ft/2V 
d,f/2 a 


distribution for the center line 


Data for the front and back surfaces of the wing were obtained from different tests 
at essentially the same stream Reynolds number. The stagnation pressure for both tests 
was 360 kN/m^, and the stagnation temperatures were 739 K and 722 K for the front and 
back, respectively. Since the value of the aerodynamic heat-transfer coefficient is not a 
strong function of the stagnation temperature, data from different tests can be used as 
input conditions provided that the initial temperatures are the same. The initial tempera- 
ture of the model was 300 K which with the adiabatic wall ^assuming an adiabatic-to-total 
temperature ratio of 0.925 obtained by using a laminar recovery factor (^Np r j and the 
Newtonian local condition on a flat plate^ and phase -change temperatures give a T of 
approximately 0.1 for these tests. 


The chordwise heat-transfer-coefficient distributions on the windward profile 
obtained by using the phase-change data reduction procedures for both the semi-infinite 
and finite slab (heat exchange at both surfaces) are given in figure 3. The distributions 
are presented in terms of the nondimensional heat-transfer -coefficient ratio h/h s , 
where h s is the theoretical heat-transfer coefficient (ref. 6, eq. (61)) for the stagnation 
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point of a 0.2-cm-radius sphere (a 0.3-m-radius sphere scaled by the same scale factor 
as the model, 0.005). 

The data obtained from these two reduction procedures differ significantly in level; 
also, the distributions diverge for the back surface. On the front surface, the procedure 
for the semi-infinite slab gives heat-transfer coefficients which exceed those obtained 
with the procedure for the finite slab by a factor of approximately 1.25 (25 percent). For 
these data the time of phase change (tp C ) divided by the thermal -interference diffusion 

time for one-half the wing thickness ^ ~ ~ varies from 6 to 16. On the back 

surface the data for the semi-infinite slab exceed the data for the finite slab by factors 
of 1.2 to 4.75 (20 to 375 percent). Here tp^t^ ^ varies from 5 to 25. However, for 
the same tp^t^ ^2 variation as the front surface (6 to 16) the data for the back surface 
of the semi-infinite slab exceed those for the finite slab by factors of 1.23 to 2.56 (23 to 
156 percent). This result indicates the influence of the time of phase change at the oppo- 
site surface and the thickness of the wing. Comparing the time of phase change and the 
thickness of wing distribution in figure 2 with the heat-transfer- coefficient distributions 
in figure 3 indicates that the difference between the results for the semi -infinite and 
finite slabs is much more sensitive to the differences in times of phase change at the 
front and back surfaces for smaller thicknesses. 

The example heat-transfer -coefficient distribution on the wing illustrates the neces- 
sity of using the numerical finite -slab data reduction procedure with the phase-change 
heat-transfer technique on midsections of thin fins and wings. However, calculating the 
actual heat-transfer-coefficient distribution along a wing profile or phase-change iso- 
therms is time consuming and expensive. Therefore, for practical engineering applica- 
tions, a graphical presentation of the results from this numerical digital computer pro- 
gram for the finite slab with heat exchange at both surfaces which will not be restricted 
to a given slab thickness or set of thermophysical properties is needed. 


Presentation of Solutions 

Identification of independent parameters . - In order to identify the correlating 
parameters and to separate the influence of some of the independent variables from the 
influence of one convective surface on another, the simpler case for the finite slab with 
heat exchange at only one surface (x = l ) is first analyzed. The closed-form solution to 
the differential equation for the one-dimensional flow of heat at the convective surface of 
this "restricted" finite slab (ref. 2) is 


T = 1 - Y - 
“ N- 


2N. 


Nu 


^Nu(^Nu + l) + 


e J 


( 2 ) 


8 



where 



a nd for j = 1, 2, . . . are the positive roots of 
V tan 77 = N Nu 

with the initial and boundary conditions of 

T(x,0) = 

dT(O.t) . 0 
dx 

dTtt.t) h [ T a» - T tQ.tj| 
dx k 


Equation (2) is plotted in figure 4 in terms of N-^ u as a function of \j/ for con- 
stant values of T. The range of the variables given should cover that encountered in any 
phase -change heat -transfer data reduction process. 

In contrast to the solution for a semi-infinite slab (appendix B, eqs. (B5) and (B7)), 
the thermophysical properties in the "restricted" finite-slab solution (eq. (2)) are no 
longer coupled as \Jpck, and time appears to the first power in the finite-slab param- 
eter i^, whereas the solution for a semi-infinite slab depends on \[t. This is not to say 
that the solution for a finite slab depends linearly on t since both solutions must agree 
for times less than the thermal -interference diffusion time of the slab. Also, the solu- 
tions for a finite slab depend on the thickness of the slab l, whereas l does not appear 
in the solution for a semi-infinite slab. 

In order to present the solutions for a finite slab in a manner more applicable to 
the phase-change heat-transfer technique, consider the independent finite-slab param- 
eters N Nu = ih A an d 1 ^= a tp C y4^. From equation (1), l is proportional to \Jat^ 
therefore, N Nu is proportional to h^t^/^pck. Comparing this parameter with the 
solution for a semi-infinite slab (eq. (2)) suggests a correlation parameter of the 
form h/h s j s (the heat -transfer coefficient for a surface of a finite slab divided by the 
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heat-transfer coefficient calculated for the same point from the solution for a semi- 
infinite slab). In the parameter \}s = , a/fi may be replaced by 0.2/t^j 

(eq. (1)); therefore, i p is directly proportional to tp^t^ thus, the parameters that 
define the restricted finite-slab solution are h/h S i S , tp C /td and T. 

A plot of the parameter h/h s j s (at surface 1) as a function of tp^t^ ^ for a 
finite slab with no heat exchange at surface 2 (eq. (2)) is given in figure 5 for three fixed 
values of T. Since h s ^ s is a function of T, the correction factor h/h s j s becomes a 
rather weak function of T. However, for the same value of the phase-change time 
parameter, the correction factor decreases (below 1 - greater departure from the value 
for a semi-infinite slab) as T increases. Also, the ratio h/h s ^ s is unity for 
tp C i/td i = 1. This result substantiates the validity of the constant in equation (1), which 
was obtained by empirical means (ref. 1). 


Effect of heating at both surfaces .- An example showing the effect of heat input 
(Nusselt number) at surface 2 on the ratio of heat-transfer coefficient for a finite slab to 
that for a semi-infinite slab at surface 1 is presented in figure 6. These curves were 
obtained from the numerical data reduction procedure for a finite slab. In figure 6 
h/hgig for surface 1 is plotted as a function of the time at which the phase change occurs 
at surface 1 nondimensionalized by the thermal -interference diffusion time of one-half 
the thickness of the slab (td,z/2 is used because heat is entering the slab at both sur- 
\ — T o ‘ Tj 

faces). Figure 6 is presented for Tj = 0.348 and — - — t = 1 and contains lines 


for constant values of the Nusselt number 


(%» - 1 ) 


^aw,l " T i 
at surface 2. 


The influence of heat 


transfer at surface 2 on the heat-transfer correction factor h/hg^g at surface 1 is 
strong for tp C i /t^ j/2 > 4 and increases substantially as this parameter increases. 

The effect of changing the front-surface Tj on the front-surface correction 
factor for a constant back-surface Nusselt number can be ascertained from figure 7. 
Here, the correction factors for a finite slab with a constant heat-transfer coefficient at 
the back surface (assuming l and k are the same for all solutions) for three values of 
T^ (t^ = 0.239, 0.348, and 0.455) are shown, along with the correction factors for a 
finite slab with no heat exchange at the back surface for the same three values of Tj. 
The finite-slab correction factor solutions for N^ u 2 > 0 depend on T^ in an inverse 
manner to that of the restricted (n^ u 2 = 0) solutions. For no heat exchange at the back 
surface (Nj^ u 2 = 0) and for a given phase-change time, the solutions for a finite slab 
move away (clockwise in fig. 7) from those of the semi-infinite slab (h/h g ^ s = l) as T^ 
increases. This is because the heat transfer to the front surface increases as Tl 
increases and the higher the heat load the greater the front-surface temperature, and 
thus the heat-transfer coefficient for this inverse method feels the effects of the finite 
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dimension of the slab. For heat exchange at the back surface (fig. 7, Nj^ u 2 = 0.139) 
and for a given phase-change time, the solution for a finite slab moves toward (counter- 
clockwise) the solution for a semi-infinite slab (h/h g ^ s = l) as Tj increases. 

At first glance there seems to be a contradiction, since the N-^ u 2=0 solutions 
are a limit of the solutions with heat exchange at the back surface (N^ u 2 = Constant) and 
figure 7 indicates that they have an opposite Tj dependence. Actually, the solutions 
with the heat exchange at the back surface (N^u 2 = Constant) are moving toward the 
solution for a finite slab (for the same T^) with no heat exchange (n^ u 2 = o) as 
increases. As Tj increases, the heat-transfer coefficient at the front surface 
increases, while for these solutions (N^ u 2 = Constant) the heat -transfer coefficient at 
the back surface remains constant. Therefore, the ratio of the heat transfer at the front 
surface to that at the back surface for the same initial driving potential is increasing 
with which is in the direction of the front- to back-surface heating ratio (°°) for the 
finite slab with no heat exchange at the back surface. Thus, as Tj increases, the two 
solutions (Nnu, 2 = 0 and N^ u 2 = Constant) converge. 

Solutions for large angles of attack .- In tunnel tests of models at large angles of 
attack, the flow on the leeward surface of a wing or horizontal stabilizer is separated, 
and the heat-transfer coefficient is less than 20 percent of that on the windward side. 

The approximate regions where I12 = hq/5 and I12 = hj are shown as bands in figure 6 
because of the wide spacing of the limited number of calculated points. For the condition 
(h2 = hi/5) and for the set of temperature parameters given in figure 6, the heat-transfer 
coefficient on the windward side is within 10 percent of the value for a semi-infinite slab 
for t pc i/t^ -j/2 = 9. However, since the I12 = hj/5 band is essentially parallel to the 
h2 '= 0 (Njju 2 = solution for the finite slab with no heat exchange at surface 2, the 
heat -transfer coefficient ratio on the windward side of the section is within 10 percent of 
the h2 = 0 solution for t pc 1/2 = 20. The similar slopes of these two lines 
(h2 = 0 and I12 = h^/5) indicate that much larger values of t pc j/t^ 1 / 2 can be attained 
with little increase in the error band. 

The question of the influence of T^ on the spread of the error band now arises. 
Correction factors for surface 1 of a finite slab are given in figure 8 for constant values 
of the Nusselt number at surface 2 as in figure 6 (T^ = 0.348) but for a T^ of 0.174. A 
comparison of these two plots shows the downward rotation of the constant Nusselt num- 
ber lines as Tl decreases (discussed previously). However, a comparison of the 
1*2 = hf/5 lines in the two plots indicates that these two lines are essentially independent 
of Tj. This, along with the fact that the I12 = 0 lines are a weak function of Tj 
(fig. 5), means that the spread in the error band (between the I12 = 0 and I12 = h^/5 
solutions for t pc ^/t^ 1/2 = 20) is also a weak function of Tj. As an example, for 
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tpc lAd 1/2 = 20 the error spread is approximately 10 percent for a Tj of 0.348 
(fig. 6) and 15 percent for a Tj of 0.174 (fig. 8). This Tj range covers that which 
will occur in most wind-tunnel tests. Also, the magnitude of the initial temperature 
potential on the leeward surface compared with that on the windward surface 
(T aw 2 “ T i)/(T aW 1 ~ TjJ will affect the error band quoted (10 to 15 percent or less) 
for tp C ij t^ i/2 = 20 even though it does not enter into the N^ u 2=0 solutions. 

For values of T aw> 2 ~ Tj smaller than that of the windward surface (fig. 9), which is 
realistic for practical test situations, the lines of N^j u 2 rotate upward toward the 
Nnu 2 = 0 solution, as a comparison of figures 8 and 9 shows. This comparison also 
shows that decreasing the initial temperature potential also rotates the h 2 = hj/5 and 
h 2 = h^ solutions upward. This upward rotation of the I 12 = h^/5 line toward the 
Nj^u 2 = 0 solution (independent of T aw 2 - TjJ decreases the error band from approxi- 
mately 15 percent (fig. 8) at tpc,l/*-d,l/2 = 20 to 10 percent (fig. 9) for a T of 0.174. 

Thus, the heat -transfer-coefficient distribution on the windward surface of the mid- 
section of wings and fins at large angles of attack can be determined within 15 percent 
for tp C i/tf} ij 2 = 20 by using the distribution of the time of the phase change and thick- 
ness of the section in conjunction with the restricted finite-slab solution (fig. 4) for the 
appropriate windward surface temperature parameter. This is restrained by the pro- 
visions that Tj be greater than or equal to 0.174 and the initial temperature potential 
for the heating rate of the leeward surface (T aw 2 - Tj) be less than or equal to that of 
the windward surface. The rate at which the error band widens as these restraints are 
violated has not yet been determined. However, the Tj restraint has been shown to be 
extremely weak, and for wind-tunnel tests the initial temperature potential, for the heating 
rate of the leeward surface T aw 2 ~ Tj is always less than that for the windward sur- 
face T aw j - T^. Therefore, there is little practical need for determining the spread of 
the error band with an increasing violation of these two restraints. 

General solutions for low angles of attack .- At low angles of attack the heat-transfer 
coefficients on opposite sides of a wing or fin are of the same order of magnitude. In this 
situation, time of phase -change distributions must be obtained on both sides of the section 
as given for the example wing profile in figure 2. Results of example calculations of the 
heat -transfer parameters for a finite slab for this more general test situation are pre- 
sented in figures 10 and 11. These figures are presented for fixed values of the corre- 
lating parameter T (t^ not necessarily equal to T 2 ) and an initial temperature poten- 
tial ratio (T aw> 2 - T^y/(T aw j - T^ of 1. This equal initial temperature potential 
restraint is generally not violated or at most violated to very little extent on wings or 
fins at zero to low angles of attack. These correlations are not restricted to a given slab 
thickness or set of thermophysical properties. 
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Figure 10 is presented for Tj = 0.45 and T 2 = 0.24. In figure 10(a), h/h s ^ s 
for surface 1 is plotted as a function of t pc l/^d,l/2- This figure contains lines for 
various constant values of the correlating parameter that specifies the time at which the 
phase change occurs at surface 2 (tp Cj 2/td,i/2)- & 1 figure 10(b), h/h S j S for surface 2 is 
plotted as a function of tpc^ftd 1/2 f° r various values of tpc,l/^d,l/2’ Figure 10(b) 
can also be used to determine the heat -transfer coefficient at surface 1 if = 0.24 
and T 2 = 0.45. Thus, these two plots (figs. 10(a) and (b)) specify the heat-transfer coef- 
ficients at surfaces 1 and 2 for Tj = 0.45 and T 2 = 0.24 provided that the thermo - 
physical properties of the slab, the thickness of the slab, and the time at which the phase 
change occurs on each surface are known. For example, with T]^ = 0.45 and T 2 = 0.24, 
•suppose a/fi = 0.2 per second, t^ ^ = 0.25 second, and the phase change on surface 1 
occurs in 4 seconds and that on surface 2 occurs in 3 seconds. Thus, tpc,l/fd,f/2 f° r 
surface 1 is 16 and tp C 2 /^d ^2 f° r surface 2 is 12. Entering these values into fig- 
ure 10(a) gives a heat-transfer -coefficient ratio h/h S jL S of 0.765 for surface 1. Again, 
tp C 2 /^d 1/2 i° r surface 2 is 12 and tp C t ^ ^2 f° r surface 1 is 16. Entering these 
values into figure 10(b) gives a heat -transfer- coefficient ratio h/h s ^ s for surface 2 
of 0.212. The large differences in the front- and back-surface correction factors 
(figs. 10(a) and (b)) are caused by the influence of the T of the opposite surface. Thus, 
for a complete solution of all test conditions likely to occur, a set of plots similar to 
those of figure 10 is needed for a range of T^ and T 2 parameters. 

Figure 11 js presented in the same format as figure 10 with = T 2 = 0.1. Since 
Ti and T 2 are equal in this figure, the (a) and (b) versions as described for figure 10 
are identical; therefore, only one plot for each set of values of T is needed in the cor- 
relations where the values of T are equal. For example, consider the y/Y = 0.45 
position on the wing profile in figure 2. The temperature parameter T is approxi- 
mately 0.1 for both sides, a/fi = 0.0863 per second, ^-d,l/2 = 0.579 second, and the 
phase change on surface 1 (front) occurred in 4.4 seconds and that on surface 2 (back) in 
6.2 seconds. Therefore, tp C w ; 2 i° r surface 1 is 7.6 and ^c.,2l^d,l/2 i° r sur- 

face 2 is 10.7. Entering these values into figure 11 gives a coefficient ratio (h/h s ^ s ^ 

of approximately 0.82. For tp C i = 4.4 seconds, h s j s is calculated from equation (B7) 
to be 61 W/m^-K, and h g is given in figure 3 as 776 W/m^-K. This gives an h/h g 
of 0.06 for the front surface, which corresponds closely to that plotted in figure 3. Again, 
tp C 2 /td 1/2 i° r sur f ace 2 (back) is 10.7 and tp C 1 / 2 is 7.6. Switching the surface 

designation subscripts in figure 8 and entering these values give an (h/hg-^g^ of 
approximately 0.5. For tp C 2 = 6.2 seconds, h S j S = 50 W/m2-K (eq. (2)), and h/h s 
is calculated to be 0.033, which matches that given for the back surface in figure 3 at 
y/Y =0.45. 
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A compilation of the results from the numerical calculations is presented in fig- 
ure 12. The correction factor h/h s i s for a finite slab is given for a and T 2 
combination range from 0.05 to 0.45 and for an initial temperature potential ratio 
(Taw, 2 " T i)/( T aw,l " T i) This range of T should cover most test conditions 

likely to be encountered in phase-change heat-transfer wind-tunnel tests. Instead of 
giving the correction factor at surfaces 1 and 2 for a T combination as in figure 10, 
these plots are given in terms of only the correction factor at surface 1. The correction 
factor at surface 1 (h/h s ^ s ) becomes larger than 1 (fig. 12) when heat is lost from sur- 
face 2. To obtain the correction factor at surface 2, one considers the plot for the 
reverse set of T, as previously stated, and treats surface 2 as if it were surface 1 in 
both the abscissa and ordinate; also, the lines of constant phase-change time at surface 2 
should be treated as if it were surface 1 . 

In using these plots, one should be aware that errors associated with small mis- 
matches in either Tp C or tp C increase as the thermal interference diffusion time 
decreases and as the difference in the time of phase change at opposite surfaces 
increases. 


Analysis of Restraints and Variable System 

In the previous section the need for determining the exact error incurred in the 
finite -slab correction factor correlation for a specific restraint violation was .stated to 
be of little value. This assessment was made on the premise that the specified bound- 
aries of the restraints are seldom violated in phase-change heat-transfer tests on thin 
surfaces. (See table I.) However, the limit within which a restraint must be contained 
in order to remain within a given error band is of interest. 

TABLE I.- BOUNDARIES OF RESTRAINTS 


Type of surface 

Solution 

Thick sections (any 0) 

Semi-infinite slab 

Thin sections (zero to 
small 0) 

Finite -slab solution with heat exchange 
on both surfaces 
Figures 9 to 13 

T combination range from 0.05 to 0.45 
(T a w,2 - Ti)/(T aw>1 - = 1 

Thin sections (moderate to 
large 0) 

Finite -slab solution with heat exchange 
at only one surface 
Figure 4 

T s 0.174 (weak restraint) 

( T aw,l “ T i) = ( T aw,2 “ ^i) 
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The heat -transfer-coefficient correction factor correlation for a finite slab 
with a phase change on both sides showing the change in the correction factor for a 
given change in the initial temperature potential ratio (T aw 2 - TjJ/(T aw ^ 

and T is given in figure 13. This figure was initially for = T 2 = 0.45 and 
( T aw,2 " T i)/(Taw,l " T i) = 1 ( case A). The initial temperature potential ratio was 
changed from 1 (case A) to 0.870 (case B) to 0.773 (case C) with the values of T 
remaining at 0.45 to illustrate the possible error in relaxing this restraint. For 
example, for tp C ^t<j i/2~ 8, the difference between h/h S j S for values of 
( T aw,2 " T i)/(Taw 1 " T i) * and 0.773 (using h/h s i s for the initial temperature 
potential ratio of 1 as the standard) is approximately 37 percent for tp C 2 /^, 1 /2 = 2 
and 24 percent for tp C 2 /^A ij 2 = 4. This type of percent deviation is shown in figure 14 
as a function of initial temperature potential ratio. The percent deviation from the 
h/hgis value for an initial temperature potential ratio of 1 decreases linearly to zero as 
the initial temperature potential ratio increases to 1 for constant values of the time of 
phase-change parameter at surfaces 1 and 2. For a given initial temperature potential 
ratio, the percent deviation at surface 1 increases substantially with an increase in the 
time -of-pha!se -change parameter at surface 1 and decreases substantially with an 
increase in the time-of-phase-change parameter at surface 2. Thus, figure 14 indicates 
that the unit initial temperature potential ratio restraint is rather strict in the sense that 
if it is violated by as much as 10 percent in certain time-of-phase-change parameter 
regimes (tp C 1/2 < &> *p C 2 /^d 1/2 < 4), the errors" which are incurred will be 
extremely large. This is modified somewhat by the fact that for wings at high angles of 
attack in hypersonic flows the lower limit of the initial temperature potential ratio is in 
the vicinity of 0.87. Also, the correlation of the type given in figure 13 should be used 
for zero to moderate angles of attack, where the initial temperature potential ratio is 
near unity. For large angles of attack, the solution for a finite slab with heat exchange 
at only one surface (fig. 4) should be used. In this situation, the error was shown to 
decrease as the initial temperature ratio decreased from unity. 

To assess the strength of the T restraint, the values of T were changed from 
0.45 (case C) to 0.59 (case E) while the initial temperature ratio was kept at 0.773. The 
effect of these changes on the correction factor is shown in figure 13. There is very 
little difference in the lines for Tj = T 2 = 0.45, (T aw ^ - T j) / (T aW) j - Tj) = 0.773 

(case C) and the lines for Tj = T 2 = 0.59, (T aw 2 - Ti)/(T aw> j - Tj) = 0.773 (case E). 
This result confirms the weakness of the T restraint for the correlations at zero to 
moderate angles of attack. 

To insure that h/h S j S is a function only of Tj, T 2 , (T aw ^ - “ Tj), 

*pc lAd 1/2’ and *pc 2 /^ 6 . 1/2’ ca l culations were made for Tj = 0.45, T 2 = 0.45, 
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and (T aw> 2 - Tj)/(T aw j - Tj] = 0.773 (case D) as in case C but with different values 
(T pc ,l ~ ^i)* (T pc ,2 ” ^i)’ (^aw,l " ^i)’ an< ^ (^aw^ " ^i) - The resu lts (identical 
curves for cases C and D in fig. 14) verify that the five variables mentioned completely 
specify the one -dimensional finite-slab system. Also, many cases which included differ- 
ences in thickness and thermophysical properties were calculated. Each correlated with 
the tpg/t^ j/2 parameter. 


CONCLUDING REMARKS 

A numerical computer program has been developed to calculate the heat -transfer 
coefficients for both sides of a one -dimensional finite slab subject to the boundary con- 
ditions ascribed to the phase-change coating heat-transfer technique. 

The data reduction procedures for a semi-infinite slab used with the phase-change 
coating technique can give large errors in heat-transfer coefficients on thin sections . 

In a typical tunnel test situation where a thin wing was exposed to heating on both sides, 
the data reduction procedures for a semi-infinite slab gave heat-transfer coefficients as 
much as 375 percent too high on the side with the lowest heating. The error at the oppo- 
site surface position, where the level of heating was approximately 5 times higher, was 
29 percent. 

The results from the procedure for a finite slab are presented in the form of cor- 
rection factors to the solution for a semi-infinite slab in terms of parameters normally 
used with the phase-change heat-transfer technique. These correlations are not 
restricted to slab thickness or thermophysical properties and are easily used to obtain 
data on thin model sections. 

At large angles of attack the heat -transfer -coefficient distribution on the windward 
surface of the midsection of a wing or fin (where lateral conduction is small) may be 
determined within 15 percent for nominal wind-tunnel test conditions by using the solution 
for a finite slab with no heat exchange on the leeward surface. This is restrained by the 
provisions that the finite-slab temperature parameter on the windward surface T^ be 
greater than or equal to 0.174 and the initial temperature potential for the heating rate 
of the leeward surface be less than or equal to that of the windward surface. The Tj 
restraint is extremely weak, and the initial temperature potential restraint is practically 
never violated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 8, 1973. 
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APPENDIX A 


COMPUTER PROGRAM FOR FINITE SLAB 

The heat -transfer computer program described in the section entitled "Analytical 
Method for Finite Slab" calculates the heat-transfer coefficients for both sides of a one- 
dimensional finite slab subject to the boundary conditions ascribed to the phase-change 
coating technique (ref. 1). 

Program Input 

FORTRAN IV NAMELIST with the name NAM1 is used to load the input data. The 
following list contains the input variables with the dimensions used in the program: 


FORTRAN 

variable 

Symbol 

Description 

TIME0 


initial time, sec 

T0 (30) 

Ti 

initial temperatures, °R 

SIGMA 

a 

Stefan-Boltzmann constant, Btu/ft^-sec-°R^ 

EMS 

e 

emissivity . 

RH0 

P m 

density, lbm/ft^ 

NKPTS 


number of entries in conductivity table (KMTAB) 

TKMTAB (10) 


temperature table for KMTAB, °R 

KMTAB (10) 

k 

conductivity table, a function of temperature (TKMTAB), 
Btu/ft-sec-°R 

NCPTS 


number of entries in specific -heat table (CPTAB) 

TCPTAB (10) 


temperature table for CPTAB, °R 

CPTAB (10) 

c 

specific-heat table, a function of temperature (TCPTAB), 
Btu/lbm-°R 
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APPENDIX A - Continued 


FORTRAN 


variable 

Symbol 


Description 

LEN 

l 

thickness of slab 



N0 


number of blocks or divisions in chord length (either 20 
or 30) 

PRFREQ 


print frequency 

T1FIN 

T pc,l 

phase-change temperature for surface 1, °R 

TNFIN 

T pc,2 

phase-change temperature for surface 2, °R 

DELTAT 

At 

delta time or computing interval, sec 

DELHI 

Ahj 

value used to increment or decrement h^, Btu/ft^-sec-°R 

DELH2 

Ah 2 

value used to increment or decrement h 2 , Btu/4t^-sec-°R 

TEST 


closeness test for values of delta temperatures, °R 

NTTPTS 


number of entries in total-temperature table (TT0TAB) 

TTTTAB (10) 


time table for TT0TAB 

TT0TAB (10) 


total -temperature table, a function of time (TTTTAB) 

TAW1RAT 


ratio of adiabatic to total temperature at surface 1 

TAW2RAT 


ratio of adiabatic to total temperature at surface 2 

RN 


nose radius of reference sphere, ft 

CPG 

c 

specific heat of test gas, Btu/lbm-°R 

PT0T 


stagnation pressure of flow 

PT2PT 


total pressure ratio across normal shock 


18 



APPENDIX A - Continued 


FORTRAN 

variable 

Symbol 

Description 

P1PT2 


ratio of static pressure to total pressure behind normal 
shock 


STIME1 


t np i time of phase change at surface 1, sec 


STIME2 


tp C> 2 time of phase change at surface 2, sec 


The only restriction on these inputs imposed by the programing procedure is that 
STIME1 must be greater than or equal to STIME2. The size of the time step (At) and the 
thickness of the slab ( 1 ) along with the block thickness distribution dictate the accuracy 
of this program for a given thermal diffusivity of the slab. Time step sizes of up to 
0.025 second were used with slab thicknesses of up to 0.015 feet to obtain heat-transfer 
coefficients accurate to within 1 percent of restricted finite-slab solutions (ref. 3) for a 
material with a thermal diffusivity of 5 x 10“® ft'Vsee. The distributions of the thickness 
of the blocks through the finite slab (built into the program) are as follows: 

If Nj2f = 20 (used for 1 < 0.005 ft with a = 5 x 10”® ft^/sec): 


Block number Thickness, w 


1,20 0.0051 

2,19 .011 

3,18 .021 

4,17 .031 

5,16 .041 

6,15 .051 

— 7,14 - .061 

8,13 .0851 

9,12 .101 

10,11 .101 


If N(2f=30 (used for 1 ^ 0.005 ft with a = 5 x 10"® ftVsec): 


Block number 

I, 30 
2,29 
3,28 

• 4,27 . : 
5,26 
6,25 

7,24 . . 

8,23 
9,22 
10,21 

II, 20 
12,19 
13,18 
14,17 
15,16 


Thickness, w 

0.0011 

.00251 

.0051 

.00751 

.011 

.01251 

.0151 

.021 

.031 

.041 

.051 

.061 

.071 

.081 

.09651 


19 



FORTRAN 

name 

BETA1 

BETA2 

H1C0MP 

H2C0MP 

TAW1 

TAWNQ 

TIME 

DELTA TIME 
H (FRONT) 

TAW (FRONT) 
H(BACK) 

TAW (BACK) 
BL0CK NJ2( 

TEMPERATURE 


APPENDIX A - Concluded 
Program Output 


Symbol 


Description 


% 


semi-infinite-slab heat -transfer parameter at 
surface 1 


^2 


semi-infinite-slab heat-transfer parameter at 
surface 2 


hsis,l 


semi-infinite-slab heat-transfer coefficient at 
surface 1, Btu/ft^-sec-°R 


hsis,2 


semi-infinite-slab heat-transfer coefficient at 
surface 2, Btu/ft^-sec-°R 


l aw,l 


aw, 2 


t 

At 


adiabatic wall temperature at surface 1, U R 


adiabatic wall temperature at surface 2, °R 


time, sec 
delta time, sec 


hi 

Taw,l 

h 2 

^aw,2 


T 


finite-slab heat-transfer coefficient at surface 1, 
Btu/ft 2 -sec-°R 

adiabatic wall temperature at surface 1, °R 

finite-slab heat-transfer coefficient at surface 2, 
Btu/ft 2 -sec-°R 

adiabatic wall temperature at surface 2, °R 

block number (refers to block distribution through 
slab, fig. 1) 

temperature of block, °R 
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PROGRAM main < I NPUT * OUTPUT * TAPE 5= INPUT *TAPE6 = OUTPUT) 

COMMON TBAR 

0 I MENS ION TO (30 ) * TEMP ( 30 ) % X ( 30 ) * TAB 1 ( 20 ) , T AS 2 ( 30 ) , B { 30 ) * C ( 30 ) < 

1 D ( 30 ) ♦ TTAWTB ( 1 0 ) * TAW1 TAB (10) ♦ T Aw2 T AB (10)# TKMT AB < 1 0 > * KMT AB ( 1 0 ) * 
2TCPTAB ( 10 ) * C PT A6 (10) ♦ THOLD 1(5)# TH0LD2 (5 ) # 

3TMUTAB ( 33 ) * MUT AB ( 33 ) « T PRN (11)# PRNO ( 1 1 ) # TT TT AB ( 1 O ) # TTOTAB (10) 

EOU I VALENCE ( TTTTAB # TTAWTR ) 

EQUIVALENCE ( NTTPTS # NT AWPT ) 

REAL LEN#KM#KMTA6 * KMT ♦ KM A V * K M 1 # KMN 
REAL MUTAB#MUS#MUW 

OATA TAB1/.005* *01 # • 02 # . 03 * • 04 * « 05 ♦ • 06 # • 085 # 4* • 1 0 ♦ *085* * 06# • 05# 

1 .04, ,03* .02* .01 * *005/ 

DATA TAB2/* 00 1**0025# .005* .0075* .01 *.01 25 *.015#. 02*. 03 #.04#. 05# 

1 .06* .07* .08* 2*. 0965* . 08 ♦ . 0 7 * • 06 * • 05 * . 04 , . 03 * . 02 * • 0 1 5 * . 0 1 25 * • 0 1 * 

• 2. 0075 *.005 *.0025, .001/ 

DATA F * G * DELTX * E 1 * E2 /O . 0 * 1 0 • * • 2 * 2* • 1 E-6/ 

DATA TMUT AS/400 • * 450 • , 500 . * 550 .*60 0. *650. * 70 0.*750 •*800**650. * 

1 90 C. *950. * 1 000 .' * l 050 • * 1 1 00 • * 1 1 50 • * 1 200 • * 1 250 • * 1 300 • * 1 350 . * 1 400 . * 

2-14 50,, 1500.* 1550. * 1 600 . * 1 650 . * 1700. * 1 750 • * 1800. *1850.* 1 900 • • 1950.* 
32000,/ 

DATA MUTAB/9 .755-6 * 1 0 .748-6 * 1 1 • 7E-6 * 1 2 • 6 E-6 * 1 3. 4 6 E-6 *14. 2SE-6* 

1 15.05E-6* 1 5.625-6* 16.57E-6* 17.29E-6* 17 .9755“ 6* 18.65E-6* 19.30E-6* 

2 1 9.925-6*21 . 1 6E~6*2l . 7SE-6 ♦ 22 . 325E-6 #22 . 86E-6 * 23 • 4 5E-6 * 23 .98E-6 * 

323. 98E-6 * 24. 50E-6 * 25. 025E-6 * 25 • 53E-6* 26 .03E-6* 26. 53E-6* 27. 02E-6 * 

4 27 • 05E-6 * 27 . 97E — 6 * 28 .4 3E — 6 * 28.885 — 6 * 29 . 33E-6 * 29. 78E-6/ 

DATA TPRn/4 00 . *45 0. * 500. * 550 . * 600 . * 650 • * 700 . * 750 • * 800. * 65 0 • * 900 •/ 
DATA PRNO/. 7 305 * • 72 1 5 * . 7 1 35 * • 706 * .7* .694 * .69* • 686* .6835 * .682 # .681 / 

-* 

INTEGER PRFREQ 

NAMELIST /NAM 1 / T I MEO * TO * S I GMA-* EMS * RHO * NKPTS * TKMT AB * KMTAB * NCPTS * 
1 J.CPT A B * CP TAB *_lEN *_N0 • PRF REQ_* „T l.F I N * J NF IN *_p_EL TAJ *j?EL H i * DEL HZ * TEST, 

3 *NTTPTS * TTTTAB .TTOTAB* TAW 1 RAT *TAW2RAT*RN* CPG.PTOT * PT2PT , P 1 P T 2 
4 , ST 1 ME 1 ♦ ST 1 ME 2 

NAMEL I ST/DEBG1 / BETA 1 *B5TA2 * HI COMP * H2CCMP * TAW1 * T A W NO 
NA MEL I ST/DEBG2/ HI *H2*TH0LD1 * TH0LD2 * TEMP 

- -NAMELIST -/DEBG3/ THOLD 1 ♦ THOLD2 * PT 1 H 1 *PT 1H2 * PT2H 1 * PT2H2-* DT 1 * D-T2 * 

1 DE T * DELHI *DELH2 

external fofb 

7C 

1 READ ( 5 * NA M 1 ) 

IF ( EOF *5)5* 10 
5 STOP 
10 CONTINUE 

WRITE (6.NAM1 ) 

I PFCT = 0 

DO 11 1=1 * NTTPT S 

T A W 1 T AB ( I ) = T AW 1 RAT * TTOTAB ( i ) 

'll T A W2T AB ( I ) = TAW2RAT* TTOTAB(l) 

IF (N0.GT.20) GO TO 20 
DO 15 • I =1 *N0 
15 X ( I ) * LEN * T AB 1 ( I ) 

GO TO 30 

20 DO 25 1=1* NO 

25 X ( I ) = LEN * TAB2 ( I ) 

30 CONTINUE 

* 

* COMPUTE HI AND H2 FROM SEM I - I NF I N 1 TE SLAB APPROX I MAT I ON 
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CALL FTLUP ( ST IME1 * TAW I t 1 n NT AWPT * TT A WTB ^TAWlTAB) 

CALL FT LUP ( S T I ME2 »TAWNOi 1 » NT AWPT *TTAWT0*TAW2TAB ) 

T6AR = ( T 1 F I N— TO ( 1 ) ) / ( T A W 1 -TO ( 1) ) 

CALL I T R2 ( BETA 1 * F * G ♦ DELTX , r OFB t £ 1 *E2*50* I CODE ) 

TEAR = ( TNF I N-TO (NO ) )/ < TAWNO-TO (NO ) ) 

CALL I TR2 (BETA2 * F DELTX * FOFB ♦ El »E2i50» I -CODE ) 

CALL FTLUP (TIFINtKMl ♦ ! * NKPTS ♦ TKMT AB t KMT AS ) 

CALL FTLUP (TNFIN,KMNil * NKPTS % TKMT AB * KMT AB > 

CALL FTLUP ( T 1 F I N , CPi ♦ 1 * NCPTS * TCPT AB * CP TAB ) 

CALL FTLUP (TNF IN* CPN * 1 *NCPTS*TCPT A8 * CPT AB ) 

H1C0MP= BET A 1 * SORT <RHO*CP 1 *KM 1 /ST I ME 1 > 

H2C0MP = 6ETA2 * SORT (RHO*CPN*KMN/ST I ME2 ) 

HI =H 1 COMP 
H2 =H2C0MP 
WRITE (6i DEBGI ) 

40 DO 130 K=lo 
DO 50 J~ 1 ♦ NO 
50 TEMP ( J ) = TO ( J ) 

TIME= TIMEO 

* GO HERE TO UPDATE TIME AND RECOMPUTE 

60 TIME=TIME+ DELTAT 
DO 70 J= 1 *NO 
B ( J)=0.0 
C ( J ) =0 . 0 
70 D ( J ) =0 • 0 

* 

# DIAGONAL ELEMENTS ARE STORED IN B ARRAY 

* UPPER DIAGONAL ELEMENTS ARE STORED IN C ARRAY 

# THE ARRAYOFLOWER DIAGONAL ELEMENTS ARE EQUIVALENT TO C ARRAY 

* CONSTANT PART OF EQ* IS STORED I D ARRAY 

* SAVE 1 CONTAINS HEAT STORED COMPUTATION 

# . SAVE. 2 CONTAINS CONDUCTION COMPUTATION 

-tf 

B ( 1 ) = —HI 

CALL FTLUP ( T I ME * TAWl * I * NTAWP T t TT A WTB * TAW 1 T AB ) 

CALL FTLUP ( T I ME * TAWNO ♦ 1 * NT A\a/PT * T T AwTB ♦ T AW2T AB ) 

D ( 1 ) = -H 1 *TAW 1 +S I GMA MS* TEMP ( 1 ) **4 
D (NO ) =-H2*TAWN0+S I GMA*EMS* TEMP (NO) **4 
CALL FTLUP ( TEMP ( 1 } , KM* 1 * NKPT S*TKMTAB*KMTAB ) 

NSTOP=NO- 1 

DO 80 I =1 * NSTOP 

CALL FTLUP ( TEMP ( 1+1 )«KMTtl * NKPTS * TKMT AB * KMT AB ) 

KMAV= (KM+KMT )*«5 
KM=KMT 

CALL FTLUP ( TEMP ( I ) * CP* 1 * NCPTS * TCP TAB * CPT AB ) 

SAVE I =RHO*CP*X ( I ) /DELTAT 
S A VE2 = 2 • 0 * K M A V / ( X < I ) +X ( I + 1 ) ) 

B ( I ) = 6(1 ) -SAVE 1 -SAVE2 
B ( I + 1 ) = - -SAVE2 
C ( I ) = SAVE2 

60 D ( I ) =D ( I ) -SAVE 1 *T£MP ( I ) 

CALL FTLUP ( TEMP (NO ) * CP* 1 »NCPTSi TCPT A3 * CPT A3 ) 

SAVE 1 = RHO*CP*X (NO ) /DELTAT 
B(NO)= B ( NO ) — S A VE 1 - H2 

D (NO ) = D ( NO ) -SAVE 1 * TEMP (NO ) 

CALL TRIDMAT ( i C*S*C*DiTfMP*NO) 

IF ( ABS (TIME-STIME2 ) ♦ GT • .00001) GO TO 85 

T HOLD 2 (K )= TEMP (NO) 
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00001 ) GO TO 60 


GO TO 60 

85 IF ( flBS (T 1ME-ST I ME 1 ) . GT . . 

T HOLD 1 (K ) =TE VIP ( 1 ) 

* 

GO TO ( 90 t l 00 i 1 1 0 i 1 £0 » 1 30 ) K 
90 HI =H 1 COMP-I-DELH 1 
GO TO 1 30 

ICO HI = H 1 COMP— DELHI 
GO TO 1 30 
110 HI =H1C0MP 

H2 =H2C0MP+DELH 2 
GO TO 130 

120 H2 =H2C0MP— 0ELH2 
130 CONTINUE 
•it- 

* T HOLD i ( 1 ) AND THOLD2 ( 1 ) CONTAIN TEMPERATURES COMPUTED AS F(H1 ,H2) 

* THOLD 1(2) AND TH0LD2<2) CONTAIN TEMPERATURES COMPUTED AS F ( H 1 + DEL H 1 , H2 ) 

* T HOL D i ( 3 ) AND THOLD 2 ( 3 ) CONTAIN TEMPERATURES COMPUTED AS FCH1-DELH1 ,H2) 

* THOLD 1(4) AND TH0LD2<4> CONTAIN TEMPERATURES COMPUTED AS F < H 1 * H2 + DELH2 ) 

* THOLDl (5)AND TH0LD2(5) CONTAIN TEMPERATURES COMPUTED AS F (HU H2 -DELH2 ) 

■>f 

PT 1 H 1 = (THOLDl (2 ) -THOLDl ( 3 ) )/ (2 .0*DELH 1 ) 

PT 1 H2= (THOLDl m ) -THOLDl (5 ) )/ ( 2 • 0*DELH2 ) 

PT2H1 = (TH0L02 (2 ) -TH0LD2 (3)1/ ( 2 • 0*DELH 1 ) 

PT2H2= (TH0L02 (4 1-TH0LD2 (5 ) ) / (2 . 0*DELH2 ) 

of 

DTI = T1FIN - THOLDl (1) 

D T 2 = TNFIN - TH0LD2 ( 1 ) 

■>c 

' OET = PT1H1* PT2H2 +PT1H2*PT2H1 

* 

DELH I- = -( D-T-1-* P -T-2 H 2 — D-T 2-*-P-T 1-H2-) /-DE-T- 

, DELH2 = (DT2*PT1H1 -QT 1 *PT2H 1 ) /DET 

IF (A5S(DTl )#LE • TEST ♦ AND* ABS < DT 2 ) *LE • TEST ) GO TO 150 
..HI COMP = H1C0MP +0ELH1 
H2C0MP = H2C0MP +DELH2 ~ 

HI' = HICOMP 
H2 = H2C0MP 
GO TO 40 

* 

1 50 DO 1 60 J=U NO 
160 TEMP ( j ) =T0 ( J ) 

time=timeo 

HI = HI COMP 
H2= H2C0MP 

170 TIME =TIME +DELTAT 
• DO 180 J= 1 < NO 
B ( J ) =0 • 0 
C ( J ) =0 • 0 
180 D ( J ) =0. 0 
B( 1 )= -HI 

CALL FTLUP (T IME-»TAW1 U • NTAWP T * TT AWT8 « T AW 1 TAB ) 

CALL FTLUP (T IME.TAWNO. I * NT AWPT • TTAWTB n TAW2TAB ) 

D < 1 ) = — H 1 *TAW1 + S I GMAtfEMS* TEMP ( 1 )**4 
D ( NO ) =-H2*TAWN0 + S I GMA*EMS* TEMP ( NO ) **4 
CALL FTLUP ( TEMP ( 1 )«KM« 1 , NKP TSiTKMTAB t KM TAB ) 

NSTOP= NO-1 
DO 190 1=1 * NSTOP 
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CALL FTLUP ( TEMP ( I + 1 ) *KMT « 1 t NKPT S • TKMTAB « KMT AB ) 

KMAV" • 5 *(KM +KMT } 

KM =KMT 

call FTLUP ( TEMP (I > ♦ CP* 1 ♦ NCPTSiT CPT AB * CPT AB ) 

SAVE 1 = RHO *CP* X ( I ) /DELT AT 

SAVEa = 2 • 0* KMAV/(X( I >+X( 1 + 1 ) ) 

B ( I ) - 6' ( r ) -SAVE 1 -SAVE2 

B ( I + 1 ) = —SAVE 2 

C ( I ) - SAVE2 

190 D(I)- D ( I ) -SAVE 1 *T£MP(I) 

CALL FTLUP ( TEMP (NO ) «CP* 1 t NCPTS * TCPT AB * CPT AB ) 

SAVE1 = RHO*CP* X ( NO ) /DELT AT 
B ( NO ) = B ( NO ) -SAVE 1 -H2 
D ( NO ) = D(NC) -SAVE 1 *TEMP ( NO ) 

CALL TR I D M A T ( C « 8 ♦ C ♦0»TtMP , NO ) 

IPFCT= IPFCT+ 1 

IF ( IPFCT.NE.PRFREQ ) GO TO 195 
I PFCT = 0 

WRITE < 6 * 1 000 )T IMS* DEL TAT * Hi * H2 * T AW 1 *TAWNO 
1300 FORMAT! 1X///6X6HTIME - F 1 0 • 4 * 6X 1 3HDEL TA TIME = F6 • 4 //7X9HH ( TOP ) = 

1 El 6.8* 5X1 2HH (BOTTOM ) = E 1 6 . 8/7X9HTA W ( TOP ) =E 1 6 • 8 * 5X 1 2HTA W ( BOTTOM ) = 

2E 16* 8//6XSHBL0CK NO 1 1X1 1 HTEMPER A TURE// ) 

WRITE (6*2000 )( I * TEMP ( I ) * I = 1 * NO ) 

2000 FORMAT (9x* I3*6X*E20. 8 > 

195 IF ( A5S (T I ME — ST IME1 ) • GT • •00001) GO TO 170 

R = 53 * 35 
G=32 • 2 

HOLD- 2.7*PT0T * PT2PT 

CALL FTLUP ( S T I ME 1 * T TO T * 1 * N T TPT S * TT T T AB * T TOT AB ) 

CALL FTLUP (TTOT * PRN ♦ 1*11 *TPRN*PRNO) 

SAVE l = .768*CPG* (PRN**-. 6 ) 

CALL FTLUP ( TTOT *MUS * I * 33 ♦ TMUT AB * MUT AB ) 

RH OS- HOLD/TTOJ 

DUEDX = SQRT ( 2 • 0*R*G * TT OT * ( 1 • 0 — P 1 PT2 ) )/RN 
5AVE2= ( (RHOS *MUS )**.4 )* SQRT(DUEDX) 

* - 

CALL FTLUP’(T1FIN*MUW* 1 * 33 * TMUT AB * MUT AB ) 

RHOW - H0LD/T1FIN 

HRS 1 = SAVE1 * SAVE2* ( ( RHOW+MUW >** • 1 ) 

CALL FTLUP (ST IME2 * TTOT* 1 * NT TPT S * TT T TAB , T TOT AB ) 

CALL r TLUP ( TTOT *PRN* 1*11 *TPRN*PRNO) 

SAVE! = .766*CPG*'(PRN**-#6 ) 

CALL F T LUP (TTOT * MUS * 1 * 33 * TMUT AB * MUT AB ) 

RHOS- HOLD/TTOT 

DUEDX = SQRT (2.0*R*G *TT OT* ( 1 * 0-P 1 PT2 ) ) /RN 
SAVE2=((RH0S *MUS )** *4 ) * SQRT(DUEDX) 

* 

CALL FTLUP ( TNF I N * MUW * 1 * 33 * TMUT AB * MUT AB ) 

RHOW- HCLD/TNFIN 

HRSNs SAVE 1 * S AVE2 * ( C RHOW*MUW > **• 1 > 

H 1 RAT = H 1 /HRS 1 
H2RAT = H2/HRSN 

-* 

WRITE (6*3000) HI *HRS1 * H 1 RAT * H2 * HRSN * H2R A T 
3000 FORMAT ( 1 X///5X3HH1 =E 1 4 .6 * 4X 1 5HH 1 ( REF • SPHERE ) =E 1 4 • 6* 4X 1 OHH 1. ( R AT I 0 ) = 
1 El 4 *6/5X3HH2=E 1 4 • 6* 4X1 5HH2 (REF •SPHERE ) =E 1 4 • 6 * 4 X 1 0HH2 < RAT I O > =E 1 4 . 6 > 
GO TO 1 
END 
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FUNCTION FOFSOETA) 

COMMON TBAR 
Y = ERF (BETA) 

C =1*0 — Y 

FOFB= ( 1 .0-EXP(3ETA**2)*C )- TBAR 

return 

END 


SUBROUTINE TRIDMAT ( A ,3 *C*D«T «N ) 

D I MENS I ON A ( 1 )tb(l ) »CU ) « D ( 1 ) • T ( 1 ) « W ( 30 ) • SV ( 30 ) ♦ G ( 30 ) 
THIS ROUTINE SOLVES A TR I D I AGONAL M ATR I X 
W ( i ) =3 ( 1 ) 

SV ( 1 ) = c ( 1 ) / B ( I ) 

G ( 1 ) = D ( 1 )/W ( 1 ) 

NM I = N — 1 
DO 100 K = 2 f N 
KM 1 = K- 1 

W(K) = B(K) - A (KM1 >*SV (KM1 ) 

I_F G0_T0_5 __ _ ^ _ 

4 SV(K)= C(K)/W(K) 

5 G ( K ) = (D(i<)“ A (KM1 )*G(<M1 ) >/W<K> 

1 00 CONTINUE 

T(N)-G(N) 

-- - DO- 200 K = T t NM 1 - • 

KK= N-K 

T ( Ki< ) = G ( KK ) - SVC KK ) *T ( KK+ 1 ) 

200 CONTINUE 
RETURN 
END 
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APPENDIX B 


SOLUTION FOR SEMI-INFINITE SLAB 


In the phase-change data reduction method (ref. 1) for the semi-infinite slab, the 
heat-transfer coefficient depends on the time required for the phase change to occur 
(tp C ), the temperature of the phase change (T pc j, the initial and adiabatic wall tempera- 
tures (Tj and T aw ), and the thermophysical properties of the model wall (\j pck ). -The 
relationship between the heat -transfer coefficient and the other parameters is deter- 
mined from the solution to the equation governing the transient one -dimensional flow of 
heat. This equation is 


8T _ 9 2 T 
at 9x 2 


(Bl) 


with the following initial and boundary conditions which most nearly describe the actual 
tunnel transient test: 


T (x, 0) = T i (B2) 

T(°°,t) = T t (B3) 


9T(0,t) _ hr T 
9x k l aw 



(B4) 


It is assumed (ref. 1) that the phase-change coating is at the surface tempera- 
ture T(0,t) and the time tp C is required when T(0,t) = T pc . Other assumptions are 
as follows (ref. 1): 

1 . The depth of heat penetration into the wall is small compared with the wall thick- 
ness and surface radius of curvature so that the wall acts like a semi-infinite slab 

(eq. (B3)). 

2. The model is isothermal before injection into the airstream (eq. (B2)). 

3. The surface experiences an instantaneous step in aerodynamic heat -transfer 
coefficient at time zero, and this coefficient is invariant with time (eq. (B4)). 

4. The thermal diffusivity a of the wall is invariant with temperature. 

The solution to equation (Bl) with the specified boundary conditions is (ref. 3) 


T = 1 - e& erfc(j3) 


(B5) 


26 


APPENDIX B - Concluded 


where 




2 r°° > 2 

erfcO) = -== \ e' A dX 

H Jp 


(B6) 

(B7) 

(B8) 


Equation (B5) is plotted in figure 15 in terms of the parameter p as a function of T. 
The heat-transfer coefficient is then determined from equation (B7). 



REFERENCES 


1. Jones, Robert A.; and Hunt, James L.: Use of Fusible Temperature Indicators for 

Obtaining Quantitative Aerodynamic Heat-Transfer Data. NASA TR R-230, 1966. 

2. Garrett, L. Bernard; and Pitts, Joan I.: A General Transient Heat-Transfer Computer 

Program for Thermally Thick Walls. NASA TM X-2058, 1970. 

3. Carslaw, H. S.; and Jaeger, J. C.: Conduction of Heat in Solids. Second ed., Oxford 

Univ. Press, Inc., 1959. 

4. Scarborough, James B.: Numerical Mathematical Analysis. Fifth ed., Johns Hopkins 

Press, 1962, pp. 199-222. 

5. Connor, L. E.; Sparks, V. W.; and Bhadsavle, A. G.: Heat Transfer Tests of the 

NASA-MSC Space Shuttle Configuration at the Langley Research Center Mach 8 
Variable Density Facility. TM 54/20-291, LMSC/HREC D162761 (Contract 
NAS 9-10506). Lockheed Missiles & Space Co., Dec. 1970. 

6. Beckwith, Ivan E.; and Cohen, Nathaniel B.: Application of Similar Solutions to Calcu- 

lation of Laminar Heat Transfer on Bodies With Yaw and Large Pressure Gradient 
in High-Speed Flow. NASA TN D-625, 1961. 


28 



Q| = h|(T aWit -T c> |) 

u Block 


no. 



q 2 = MW "Tc,2^ 


Figure 1.- Geometry of finite slab. 




Chordwise profile of wing on straight-wing orbiter at 69.3 percent 
Ln with time of phase-change and thermal diffusion distribution. 
i2 cm. 




Figure 3.- Result from finite- and semi-infinite-slab heat-transfer coefficient 
procedures for wing on straight -wing orbiter at 69.3 percent semispan. 

Y = 2.042 cm. 




Figure 4.- Heat-transfer solution at convecting surface of finite slab with no 
heat exchange at back surface (plot of eq. (2)). 
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Figure 6.- Front-surface correction factor for finite slab with heat exchange at 
back surface. Tj = 0.348; (T aw> 2 _ r ^'i)/( r ^'aw,l “ ^i) = ^ 






heat exchange at back surface. 
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(a) Surface 1. 

Figure 10.- Heat-transfer-coefficient correction factor for finite slab with phase change on both sides. 

T, = 0.45; T 2 =0.24; (T aw>2 - - Tj = 1. 
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Figure 11.- Heat -transfer-coefficient correction factor for finite slab with phase change on 
both sides. Tj = T 2 = 0.1; (T aw 2 - Tj)/(T aw ^ - Tj) = 1. 
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